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Abstract 

■We report ab initio calculations of the X and B -E^ states of the Bcj dimer Full valence configuration interaction 
calculations were performed using the aug-cc-pVnZ basis sets and the results were extrapolated to the CBS limit. Core-core, 
core-valence efi'ects are included at the CCSDT/MTsmall level of theory. Two local minima, separated by a large barrier, 
are found in the expected repulsive B state. Spectroscopic constants have been calculated and good agreement is found 
with the recent measurements of Merritt etal. Bound vibrational levels, transition moments and lifetimes have also been calculated. 
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Introduction 

In recent years, there has been a lot of interest in ultracold 
atom-ion scattering Ol in the atomic and molecular physics 
community. The experimental realization of Bose-Einstein 
condensation (BEC) has led to numerous applications involv- 
ing charged atomic and molecular species. The cooling and 
'trapping of such charged gases at sub-kelvin temperatures 
is a topic of growing interest. The phenomena of charge 
transport like resonant charge transfer yD and charge mobihty 
fll at ultracold temperatures have also been studied in detail. 
Other emerging fields of interest include ultracold plasmas jstl, 
ultracold Rydberg gases and systems involving ions in a 
BEC 011. 

It is well known that the Be2 dimer is a difficult problem for 
computational quantum chemistry, due to the 2s - 2p near 
degeneracy in the beryllium atom|9]. This near-degeneracy 
problem also arises in the less studied Bcj dimer. The ground 
X state is well described by a Icr^lcrllo-^lo-^^ reference, 
but the B state has a multi-reference character, as previ- 
pusly discussed by Fischer et al. jlOll . In this paper we present 
calculations on the Be^ dimer that should be useful as a starting 
point for further studies in ultracold atomic and molecular 
physics. 

We begin by describing the methods used in our calculations 
followed by a discussion of the results which include the 
potential curves, spectroscopic constants, dipole moments, 
hfetimes of the bound vibrational levels of the ^"Lt state and the 
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Figure 1: [COLOR ONLINE] Dominant molecular orbital configurations in 
the B ^1+ state. The curves in black shows the SCF curve crossings due to the 
change of the occupied molecular orbitals, whereas the curve in red shows the 
calculated FCI potential curve for the B -Z+ state. 



analysis of long range behavior and determination of the Van 
der Waals coefficients. The B state of Be^ has a double 
minima instead of a purely repulsive nature as one would 
expect. We have characterized both the inner (deep) and outer 
(shallow) well. 



Methods 
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As mentioned above, the B state cannot be adequately 
described by a single molecular orbital configuration. At 
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short internuclear separation, the dominant configuration 
is lcr^lcr^2cr^3cr^ while at large separation it becomes 
lcr^lcr^2cr^2cr^. This behavior is shown in Fig. 111. Prelimi- 
nary calculations done at the CCSD(T) level of theory find a 
kink in the potential curve for the B state at the SCF curve 
crossing (see Fig. 121). Valence full configuration interaction 
(FCI) calculations were found to give a smooth potential curve. 

Therefore, our computational approach is to perform valence 
FCI using the augmented correlation consistent polarized va- 
lence n-tuple zeta (aug-cc-pVnZ) basis set of Dunning. Ml ill We 
then extrapolated the results from the aug-cc-pVQZ and aug- 
cc-pV5Z calculations to the complete basis set (CBS) limit. 
We have used Schwenke's linear formula lll2ll to extrapolate the 
SCF energies. For extrapolating the FCI correlation energies we 
have used the following formula given by Helgaker 11311 : 



(1) 



where X, Y are 4, 5 corresponding to the aug-cc-pVQZ and 
aug-cc-pV5Z basis sets. The total valence energy is the sum 
of the extrapolated SCF and full CI correlation energies. 
Core-core (CC) and core-valence (CV) correlations were 
calculated as the difference between all-electron and frozen 
core CCSDT Q calculations done with Martin's MTsmall 
basis set The MTsmall basis set consists of a completely 
uncontracted cc-pVTZ basis set augmented with two tight 
d and one tight f functions. The calculated potential energy 
curves are corected for the effects of basis set superposition 
error by the counterpoise method of Boys and Bernardi| 16]. 
The CBS extrapolation increased the well depths of the X 
and B states by ~ 40 cm"', however Do for the outer 
well in the B state was unchanged. Scalar relativistic 
corrections were estimated to be ~ 10 cm"' and are neglected. 



The FCI calculations were done using the MOLPRD 2009 . 1 IVI 
and PSI3 01811 electronic structure programs running on a Linux 
workstation (2 quad core Intel Xeon E5520 CPU). The core- 
core and core-valence corrections were done with the multi- 
reference coupled cluster (MRCC) program lOofM. Kallay. 
Le Roy's LEVEL program f20] has been used to calculate the 
bound vibrational levels, Frank-Condon factors and Einstein A 
coefficients. Using these Einstein A coefficients we were able to 
calculate the lifetimes of all vibrational levels of the ^l,t state. 



Results and Discussions 

Potential Curves and Spectroscopic Constants 

Fig. ||3] shows the ab initio potential curves for the lowest ^E^^ 
and states of Bej dimer. We have used a standard Dunham 
analysis ll2lll to calculate the spectroscopic constants (Table [TJ. 




Figure 2: [COLOR ONLINE] The curves in black sliow a CCSD(T) calcula- 
tion done with aug-cc-pV5Z basis set, whereas the curves in red are a full CI 
calculation using the same basis set. The inset shows a discontinuity (black 
line) in the CCSD(T) curve for the B ^Zj state due to the change in reference 
configuration at larger internuclear separation. 



The bond dissociation energy (Do) is 16172 cm"' is in very 
good agreement with the experimental data of the BeJ dimer 
[22]. The X state supports approximately 70 bound 
vibrational levels. 

The B state which was expected to be repulsive has two 
minima instead. Both these wells support bound vibrational 
states. The outer well has 12 bound levels which are long Uved 
(~ ms). 



Table 1 : Calculated spectroscopic constants of Be^ 



State r, (A) B, (cm ) ai, (cm ) difXe (cm ') Do (cm ') 



X^S* 2.223 0.756 
Exp. [221 

B^Y-liln) 2.123 0.829 

B^S+ (Out) 7.106 0.074 



525.299 
498(20) 
547.452 
33.703 



4.454 

11.681 
3.548 



16172 
16072(40) 
2550 
69 



' The aug-cc-pV5Z basis was created by adding diffuse primitives with the 
following exponents to the published cc-pV5Z basis: s 0.013777, p 0.007668, 
d 0.0772, f 0.01375, g 0.174, h 0.225. 
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Figure 3: [COLOR ONLINE] The figure shows the calculated ah initio poten- 
tial curves of BeJ . The inset is a magnification of the shallow long-range well 
in the B ^Sj state (in red). The positions of the first four bound vibrational 
levels are shown (in blue) for both states of Be^ dimer Note that the energy 
scale for the inset is in cm"' . 




3 6 9 12 15 18 

R (a.u.) 

Figure 4: [COLOR ONLINE] The figure shows a plot of the computed elec- 
tronic dipole transition moment /.(xb coupling the B to the X ^Zj state. The 
dotted line (in blue) shows the a plot of R/2. 

in the inner well (~ 10"^ s). Our results for v' = - 3 agree well 
with the results of Fischer et al. ifioll . 

Table 2: Radiative lifetimes of the vibrational levels of the B ^I,* state (in s). 



Transition Moments and Lifetimes 

To compute the transition moments coupling the X and B 
states of Bcj dimer, we have used a 16 orbital complete 
active space self consistent field (CASSCF) wavefunction as a 
reference for performing multi-reference configuration interac- 
tion (MRCl) calculations. The transition moment for electric 
dipole transitions is defined as, 

fiXB(R)^{B\er\X), (2) 

where \X) and \B} are the electronic wave functions corre- 
sponding to the states X and B '^I,'^ when the two Be cores 
are separated by the distance R. Fig. H) shows a plot of the 
computed electronic dipole transition moment between the B 
^'Zg and the X ^E^^ ground states of Bcj . The transition moment 
fixB asymptotically follows the classical dipole behavior, 
fixB ~ R/2. The curve shows a zero-crossing at around 5.5 
bohrs which is approximately the same distance at which the 
dominant molecular orbital configuration changes from (Icr^ 
lo-l 1(t\ 3o-g) to (10-2 y^i 2cr^ 2(tI) in the B state of BeJ 
(Fig. ID). 

The calculated potential curves and the electronic transition 
dipole moments were used as input to Le Roy's LEVEL pro- 
gram to calculate the Einstein A coefficients coupling the vi- 
brational bound levels of the B state to the X ^2+ state. We 
have also calculated the radiative lifetimes (Table: ||2l) of the 
vibrational levels in the B state using these Einstein A co- 
efficients. Note ffiat the bound levels in ffie shallow outer well 
are extremely long-lived (~ 10"^ s) in comparison to the levels 
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Long Range Coefficients 

For large internuclear separations, the standard long-range form 
of the intermolecular potential is: 

v v («2/2 + Ce) 

Vlr{R) = Vco Y4 ^6 - ' 



Table 3: Long Range Coefficients for both X^ZJ and B^'L't states (in a.u.) 
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where Egxch is the exchange energy contribution and Voo is 
the energy of the atomic asymptote (which we have set to 
zero). Note that is the static dipole polarizability, q'2 is the 
quadrupole polarizability and C(, is the dispersion coefficient. 
Eexch is repulsive (plus sign in Eq. 131) for the state and 
attractive (minus sign in Eq. 131) for the state. 

All the parameters in Eq. 13 are common for both the X ^1,'^ 
and B ^I.'^ states. Neglecting higher order terms in Eq. Q, if 
we add the potentials for both states, the exchange term cancels 
and we get, 



(«i/2) + 



(a2/2 + Ce) 



(4) 



Table 13 lists the values of the long-range coefficients that we 
obtained from fitting our long-range data to Eq.Q. We get ai 
= 38.12 a.u. which is in excellent agreement with previous re- 
sults ll23ll . We have also calculated the quadrupole polarizabil- 
ity from a finite-field calculation using MOLPRD which gives 02 
= 300.01 a.u. Thus from the fit we were able to extract the 
value of the dispersion coefficient C(, = 124.22. This is in good 
agreement with unpublished results of Mitroy ll2411 . 



Concluding Remarks 



Accurate ab initio calculations have been performed on the X 
and B ^EJ states of the Bcj dimer Since the state has a 
shallow well near 13.4 bohr, it was necessary to include diffiise 
functions in the basis sets to describe the well accurately. Large 
augmented basis sets of the Dunning correlation consistent 
series were thus chosen and the results were also extrapolated 
to the complete basis set limit. We have corrected our valence 
only FCl results for core-core and core-valence effects using 
CCSDT calculations with both full and frozen core using 
Martin's MTsmall basis set. 
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Since the state has not been experimentally observed we 
were unable to compare our theoretical values for dissociation 
energies or spectroscopic constants [22]. However there are 
recent experimental results for the bond dissociation energy 
and (L)g for the state which are respectively 16072 + 40 
cm"' and 498 + 20 cm '. These values compare well with our 
calculated results of 16172 cm"' and 525 cm"' . 
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